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SOLUTION-GENERATING TECHNIQUES 



For stationary axisymmetric space-times many solution-generating tech- 
niques are known. These are associated with the two KiUing vectors, and 
with the internal symmetries of the Ernst equation. Colliding plane wave 
solutions also have two Killing vectors, and the main field equations can 
be transformed exactly to the Ernst equation. The main difference is that 
the solutions of the Ernst equation now contain the metric functions ex- 
plicitly rather than the field potentials as in the stationary axisymmetric 
case. It is therefore to be expected that many of the familiar solution 
generating techniques can be used to obtain new colliding plane wave 
solutions from already known, or 'seed', solutions. 

In this chapter we will consider only vacuum solutions. This will 
enable us to continue to concentrate on the collision of gravitational waves. 
Various techniques will be discussed here, and the main solutions that 
have been obtained using them will be described in the next chapter. 

Many techniques are also known by which solutions of the Einstein- 
Maxwell equations, and other non-vacuum solutions, can be generated 
from known vacuum solutions. These will be discussed later in Chap- 
ter 15. 

12.1 The colinear case 

Consider first the case when the approaching waves have colinear polar- 
ization. The metric can be taken in the form of the line element (6.20) 
with W = 0. In this case, given any solution Uq-, Vq, Mq of the field 
equations (6.22), then another solution of the same equations again with 
W = is given hy U,V and M , where 

U^Uo^-log{f{u)+g{v)) 

V = Vo + aUo (12.1) 
M = Mo + aVo+ ^^a^Uo 

where a is an arbitrary constant. 

This transformation is well known (Halilsoy, 1985). It can clearly 
be seen that, if the initial solution satisfies the boundary conditions ap- 
propriate for colliding waves, then so does the new solution. It can also 
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be seen that the transformation (12.1) is equivalent to adding to V the 
Stoyanov solution (10.29). 

It may also be observed that the transformation (12.1) can be used to 
obtain the soliton solution of Ferrari and Ibanez (19876) described in Sec- 
tion 10.4, using as seed the Khan-Penrose solution. Ferrari and Ibaiiez, 
however, originally obtained their solution using the soliton techniques 
mentioned in Section 12.6. 

Another transformation which maps one colinear solution of (6.22) 
on another, whether or not they satisfy the boundary conditions, is given 

by 

V = bVo 

M = {b'^- 1) log/V + ^(^^ - + b^Mo + const "^^^"^^ 

where b is an arbitrary constant and — f{u) + g{v). 

When b ^ ±1, it can be seen that the transformation (12.2) is singu- 
lar on the boundaries of the interaction region where f = and ^r' = 0. 
It follows that the necessary boundary conditions cannot be satisfied for 
both solutions, if the same expressions for / and g are used. If the seed 
solution satisfies the boundary conditions as described in Chapter 7, with 
Uo — — log(/o + go) and Mq containing the terms 

Mo = . . . + /coi log(| - /o) + ko2 log(i -go) + ..- (12.3) 

with 

/o = 1 - (Col«)"°^ + . . . , go=l- (Co2^)"°^ + . . . (12.4) 

then the new solution containing the terms (7.10) will also satisfy the 
boundary conditions if / and g are replaced by 

/ = I - icur^ + . . . , g=^- {C2vr^ + ... (12.5) 

where 

1 - — = 62 fl - — V i-L^b^(i-±.] (12.6) 

provided ni > 2 and n2 > 2. These equations give the new expression 
for U in terms of the null coordinates u and v. It may be noted that the 
inequalities (12.6) provide strong constraints on the range of permissible 
values of the parameter 6 in (12.2). 

It can thus be seen that this transformation with b'^ ^ 1 may be used 
to change the profile of the approaching waves. For example, it could 
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be used to derive the Szekeres solution described in Chapter 9 from the 
Khan-Penrose solution. The two transformations (12.1) and (12.2) may 
also be used in conjunction. 

It may also be noted that the transformations for V in (12.1) and 
(12.2) may be restated in the form that, if a real is a solution of Ernst's 
equation (11.18), then 

z ^{i-t^y/'^ii- z^Y'^Zo (12.7) 

and 

Z = Z^ (12.8) 

are also real solutions with arbitrary constants a and 6, though as ex- 
plained above, the possible values of h are constrained by the boundary 
conditions. 

These results are, of course, implicitly contained in the general ap- 
proach to colinear solutions described in Section 10.1. The point to note 
is that, when the approaching waves have constant aligned polarization 
so that Z and E are real and = 0, the main field equations can be 
linearized. In this general superposition of solutions is possible, as 

illustrated in (10.16) together with (10.20). Further generating techniques 
for this case have been described by Kitchingham (1984). 

In the more general case, however, Z and E are complex, and the 
field equations are essentially non-linear. In this situation, the generating 
techniques described in the remainder of this chapter are particularly 
significant. 

12.2 Rotations and Ehlers transformations 

Consider first the well known result that, if Eq is a solution of Ernst's 
equation (11.14), then 

E = e'"Eo (12.9) 

where o; is a constant, is also a solution. 

In terms of the alternative function Z the equivalent result is that, 
if Zq is a solution of (11.8), then another solution is given by 

cos § Zo — i sin § 
COS f — z sin ^ Zo 

By applying the coordinate rotation 

a; + zy = e^"/2(a;o + zyo) (12.11) 



(12.10) 
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to the line element (11.12), it may be shown that this rotation is exactly 
equivalent to the transformation (12.9), indicating that this transforma- 
tion may be simply interpreted as a global rotation of coordinates. 

In terms of the metric functions of the line element (6.20), this rota- 
tion can be shown to be equivalent to the result that, if Uq, Vq, Wq and 
Mo are solutions of (6.22), then another solution is given by 



U = Uo 

tanh V = cos a tanh 14 -|- sin o; sechl4 tanh Wq 
sinh W = cos a sinh Wq — sin a sinh Vq cosh Wq 
M = Mo. 



(12.12) 



It may be noted that, in (12.12), the function M is unaltered. Thus 
the boundary conditions are not affected, and so neither are the functions 
f{u) and g{v). In fact, it is clear that (12.12) does not generate a new 
solution, but merely a rotation of the original. 

It can also be shown that, if the seed solution is colinear so that we 
can put Wo = 0, then the transformation (12.12) reduces to the transfor- 
mation suggested by Ray (1980) as corrected by Halilsoy (1981). 

Another well known result is that, if Zo is a solution of Ernst's equa- 
tion (11.8), then 

where c is a real constant, is also a solution. This is loosely referred to as 
an Ehlers transformation (Ehlers, 1957). 

Equivalently, if Eo is a solution of (11.14), then another solution is 
given by 

l + E= ^^^^ . (12.14) 

l-f(l + i?o) ^ ^ 

In terms of the metric functions in (6.20), we may start with the 
solutions Uo: Vo, Wo and Mo of (6.22). After integrating (11.21), it can 
then be shown that another solution is given by 

U^Uo 

g2V ^ g2Vo ^ 2ce^° tanh Wo + 

(12.15) 

sinh W = sinh Wq + c cosh WqC ° 
M = Mo. 

Since the metric function M is unchanged, it follows that the bound- 
ary conditions for the new solution are satisfied if, and only if, they are 
satisfied for the seed solution. 
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The above transformation is also related to the well known result 
that, if Zo is a solution of Ernst's equation (11.8), then another solution 
is given by 

Z = Zo + ih (12.16) 

where 6 is a real constant. Equivalently, if Eq is a solution of (11.14), 
then another solution is given by 



l + f(l-K) 



l-E= ^ , ^ (12.17) 



In terms of the metric functions in (6.20), this leads to the transfor- 
mation that, if Uo, K>) Wo and Mq are solutions of (6.22), then another 
solution is given by 



U = Uo 

^-2V ^ ^-2Vo ^ 26e-^° tanh Wo + h"" 
sinh W — sinh Wq + h cosh WqC^" 

M=Mr, 



(12.18) 



This can clearly be seen to be identical to (12.15) but with the sign of V 
interchanged. 

All of the above transformations can be shown to be obtainable from 
a general rotation and rescaling of the coordinates in the x,y plane. This 
can be described by subjecting the Killing vectors and dy to any 
SL(2,i?) transformation 



dx\ ( CL h\ ( dj; 



dy J \c d j \dy 



where ad-hc=l. (12.19) 



This leaves U and M invariant but, by considering the line element (11.6), 
the Ernst function Z can be shown to transform as 

^^^(a|^_ (12.20) 

The transformation (12.20) clearly incorporates (12.10), (12.13) and 
(12.16) which include both pure rotations and those referred to as Ehlers 
transformations. When applied to metric functions on its own this trans- 
formation clearly has little physical significance. It follows that these 
transformations on their own can not be used to generate genuinely non- 
colinear solutions from colinear ones. However, they may be combined 
with some other method to obtain such solutions. 
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In particular these transformations may be applied to the potentials 
for the metric functions that will be described in Section 12.4. Strictly 
speaking it is only these transformations, applied in the potential space, 
that should be referred to as Ehlers transformations (Ehlers, 1957). It is 
then possible that a sequence both of coordinate transformations (12.19) 
and Ehlers transformations (12.20) in the potential space can be applied 
alternately to generate a further class of solutions with an infinite number 
of parameters. 

12.3 Geroch transformations 

Geroch (1971) has presented a general technique for generating solutions 
of Einstein's source-free equations from known solutions which possess a 
Killing vector As originally presented, the method is to consider a 
vacuum space-time with metric g/^i,, and to solve the equations 

for G, and subject to the following conditions: 

e^^. = F, eP^. = G, eQ^. = F^ + G^-l. (12.22) 

(A solution of these equations is known to exist.) Then, the basic result is 
that a new solution of the source-free equations, depending on an arbitrary 
real parameter a, is given by 

g'^, = FF-' {g^, - + Frj^rj^. (12.23) 

where 

F = F ((cos a - G sin a)^ -t- sin^ a) ""^ 

^ ^ (12.24) 

— F + 2 sin a cos aP^ — sin^ aQ^. 

It can be seen that F and G are respectively the norm and twist potential 
of the Killing field 

This technique determines a new solution with one extra parame- 
ter a. A second application of the method simply yields another member 
of the same one-parameter family of solutions. It may be observed that 
the new space-times also have the same Killing vector as the seed solution. 



(12.21) 
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It may also be shown that these transformations include both the 
Ehlers transformation and the rotation described in previous sections.^ 

The case when there exist two commuting Killing vectors, which 
is the case for colliding plane waves, has also been considered (Geroch, 
1972). Clearly, the above technique may be applied using any linear com- 
bination of the Killing vectors. Successive applications of the method 
using different combinations of the Killing vectors yield metrics which 
depend on the order of the vectors used. The resulting class of solutions 
depends on two arbitrary functions, and retains the same pair of com- 
muting Killing vectors. The algebraic structure of these transformations 
has been analysed by Geroch (1972). 

The Geroch transformations, which relate one vacuum solution to 
another, in fact form an infinite-dimensional group. The first really useful 
realization of the Lie algebra of the Geroch group was formulated by 
Kinnersley and Chitre (1977, 1978a,6), who demonstrated the action of 
the infinitesimal elements of the group in terms of an infinite hierarchy of 
potentials. Kinnersley and Chitre, and also Hoenselaers, Kinnersley and 
Xanthopoulos (1979), exploited this approach to derive new stationary 
axisymmetric solutions. This formalism has been extended to space-times 
with two space-like Killing vectors by Kitchingham (1984), although he 
initially applied it only in the cosmo logical context. 

Hauser and Ernst (1979, 1980) introduced a realization of the finite 
elements of the Geroch group. They also showed that the Kinnersley- 
Chitre transformations can be carried out by solving an appropriate ho- 
mogeneous Hilbert problem. 

A simpler realization Go (Si) of the Geroch group of transformations 
which is specifically applied to colliding wave solutions has been found 
by Ernst, Garcia-Diaz and Hauser (19876). With this, colliding wave 
solutions with any desired number of parameters may be constructed. 

When applying the original transformation (12.23) to colliding plane 
waves, it is convenient to start first with the Szckcrcs line element (6.20), 
in which the two Killing vectors are ^Wf^ = and = s^. It is then 
appropriate to consider a general Killing vector 

= aS^ + b6^ (12.25) 

where a and b are arbitrary constants. It follows immediately that 

F = e-^ {a^e^ coshW - 2absmhW + b^e-^ coshW). (12.26) 



For a general introduction to Geroch and related transformations see Xan- 
thopoulos (1985). 
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It is then possible to put two of the components of and each equal 
to zero. Putting 

= (0, 0, P2, Ps) , = (0, 0, Q2, Qs), (12.27) 

the remaining components of P^ must then satisfy the equations 

P2,„ = -e-^ {ae^iWu - K sinh cosh W^) - b{Uu + Kcosh^ W)) 

P^ ^ = {ae^iWy - Vy sinh cosh W^) - b{Uy + K cosh^ W)) 

P^.^ = -e"^ {aiUu - Kcosh^ W)-be-^(Wu + VuSmhW coshW)) 

Pg „ = e"^ {a{U^ - K cosh^ W) - he'^ (Wy + K sinh W cosh W^)) . 

(12.28) 

These equations are automatically integrable in view of (6.22). It is then 
possible to put 

G = aP2 + bPs (12.29) 
and to choose the components of such that 

aQ2 + bQ3^F^ -1. (12.30) 

Equations (12.21) are then satisfied. 

The new metric is then given by (12.23). The expressions for the 
new functions are rather involved. However, it is clear that 

e-^ = ((cos a-Gsinaf + sin^ a) e"^" (12.31) 

which indicates that the boundary conditions are satisfied for the new 
solution, if they were satisfied for the seed solution. 

In the case when the initial solution has colinear polarization, this 
method can be considerably simplified. Such a situation was considered 
by Panov (1979&), and will be discussed further in Section 13.2. 

12.4 The Neugebauer— Kramer involution 

Chandrasekhar and Ferrari (1984) have shown that the main field equa- 
tions for colliding plane waves can be written in the form of Ernst's equa- 
tion (11.18) involving the complex function Z. In this case, the real and 
imaginary parts of this function, denoted by Z = x + ^ire the metric 
coefficients given by (11.2-4). Using these components, the main field 
equations can be written as 
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and 

^-.).-(^^-.),. = 0. (12.33) 

These are respectively the real and imaginary parts of (11.18). 

The second of these equations immediately implies that there exists 
a function ^> such that 

= ^-^^^t, = (12.34) 

It follows that $ must satisfy the equation 



At this point, it is convenient to introduce a new function defined 

by 

= Vl-tVl-'Z^ X"^- (12.36) 

From (11.1) and (10.11), it may be seen that —'^ is the coefficient of dx^ 
in the line element (11.4). 

With this definition and (12.34), the equations (12.32) and (12.35) 
can now be rewritten in the form 



and 



*2 ^^J-^-^^^^^J- = 0- (12-38) 

These equations may be seen to be identical to (12.32) and (12.33) 
with \I/ and $ replacing x and lo. Thus, ii = x -\- ilo '\% solution of 
(11.18), then another solution is given by 

Z = * + z$ (12.39) 

where is given by (12.36), and $ is chosen to satisfy (12.34). This 
correspondence is known as the Neugebauer-Kramer involution. It is 
known to map real metrics to real metrics only for space-times with two 
space-like Killing vectors. In modern terminology, it is clear that the 
equations (12.34) provide an auto-Backlund transformation for the Ernst 
equation (11.18). 
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In order to obtain the remaining metric function M it is necessary 
to integrate equations (7.9), which may be rewritten in the form 

,/ ^ (12.40) 

which is equivalent to (11.21). In terms of the alternative functions \E' 
and ^» these become 

^ 2(/ + o) * 2 *2 ^^/^^/^ 

^■1 (12.41) 

It follows from this that, for any solution with = x + io; and Mq, 
there exists another solution with Z = + i$ = Zq and M where 

M = Mo + ^C/o+logx. (12.42) 



This immediately indicates that, if the initial solution satisfies the bound- 
ary conditions for colliding plane waves, then so also does the new solu- 
tion, provided 

U = Uo. (12.43) 

A second application of this technique only leads back to the original 
solution, apart from an arbitrary complex constant. The presence of 
such a constant has already been considered in Section 12.2 following 
equation (12.16). 

It may be noted in passing that, in the colinear case in which a; = 
and $ = 0, this generation technique is equivalent to putting V = Uo — Vq 
which is effectively contained in (12.1) but with a change in the sign of V . 
This, for example, would generate the Ferrari-Ibahez solution from the 
Khan-Penrose solution. 

It is of interest to observe that, for stationary axisymmetric space- 
times, the main field equations can be written in the form of Ernst's 
equation (11.8), whose solution Z \E' -|- i$ is considered as a complex 
potential for the metric functions x c^nd u that are obtained from (12.36) 
and (12.34). In contrast, for colliding plane waves, both the potential 
(12.39) and the metric function (11.5) satisfy the same equation. Thus, in 
this case, the Neugebauer-Kramer involution is equivalent to a generation 
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technique in which any solution has a dual solution in which the metric 
functions and the associated potentials are interchanged. 

12.5 A combined transformation 

In a recent paper, Halilsoy (1988fc) has considered a technique which he 
applied initially in the context of colliding shock electromagnetic waves 
(Halilsoy, 1988a). This can be described as follows. 

If Vo is a solution of the vacuum equations (9.3), (10.2) or (10.13) 
corresponding to a collision of colinear gravitational waves (with Wq = 0), 
then a noncolinear solution is given in terms of the function E satisfying 
(11.14) by 

E = Pe'^ (12.44) 

where 

p2 cosh aVo — cos a 

cosh aVo + cos a (12.45) 
tan $ = — tan a coth oVq 



where a is an arbitrary constant. The original solution is obtained when 
a = and a = 1. 

In terms of the alternative function Z, the original solution may be 
described by a real function Zq and the new solution satisfying (11.8) is 
then given by 

ztanf + Z« 

Z = ^ ^ ° . 12.46 

l-ztanfZ« ^ ^ 

This can thus be seen to be a transformation of the type (12.2) followed 
by an Ehlers transformation in a slightly more general form than (12.13). 

In terms of the metric functions of the Szekeres line element (6.20), 
the original solution is described by the three functions Uo, Vo and Mq 
which are each functions of /o and go where, as in (7.8), it is convenient 
to put 

e-^° =fo + 9o, e-^° = fl^ ^-s. (12.47) 

V Jo + Qo 

and, as in (12.4), 

/o = ^ - (Colw)"°^ + . . . , go=l- (Co2^)"°^ + . . . . (12.48) 
The new solution (12.44) is given by 



tanh V = cos a tanh aVo 
sinh W = tan a cosh aVo- 



(12.49) 
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With this, equations (7.7) can be integrated to give 

M = (o^ - 1) log/X + - ^)Uo + a^Mo + const (12.50) 

which is similar to (12.2). Now, using the same arguments as in Sec- 
tion 12.1, the new solution will only satisfy the appropriate boundary 
conditions if the forms of the functions /o and go are modified. Accord- 
ingly we now put 

e-^ = f + g, e-^=-^e-^. (12.51) 

Vf + 9 

The parameters / and g are now expressed in terms of the null coordinates 
u and V in the form 

/ = i - {c^ur^ + . . . , ^ = i - {C2vr- + ... (12.52) 

where 



(12.53) 



With this restriction, the boundary conditions are satisfied, and 

S = a^So. (12.54) 

If, in the original solution, the functions /o, go and Mq have been 
chosen such that the terms fog'o in (12.47) have already cancelled the 
terms koi log(| — /o) and ko2 log(| ~ 5*0) in So, then the new solution for 
M will be given by 

M = |(a2 - l)Uo + a^Mo. (12.55) 

The method described in this section has been used by Halilsoy 
(19886) to obtain a class of solutions which appears to generalize the 
Szekeres solutions. However, since it has been obtained by a real trans- 
formation followed by an Ehlers transformation, as described at the end of 
Section 12.2 it is essentially simply a rotation and rescaling of the Szek- 
eres solutions. It does not therefore describe the collision of genuinely 
non-aligned gravitational waves. 



12.6 Other methods 



Within the context of stationary axisymmetric space-times, many tech- 
niques have been developed in recent years by which new solutions may 
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be obtained from initial vacuum solutions.^ Most of these methods can 
be apphed also to the colliding plane wave situation. The main difference 
being that, in this case, the solutions of Ernst's equation may themselves 
either be metric functions or potentials for those functions as described 
in Section 12.4. 

The other difference relates to the change of boundary conditions. 
For stationary axisymmetric solutions, it is usual to impose the condition 
that the space-times are asymptotically flat, and techniques for generating 
new solutions of Ernst's equation in which the space-time satisfies this 
condition have been developed. One of the most general of these is that 
of Cosgrove (1977), which generalizes the Tomimatsu-Sato solutions to 
arbitrary continuous deformation parameter 5. However, for colliding 
plane waves, the boundary conditions are of a totally different character, 
and many of the solutions of Ernst's equation that are appropriate for 
axisymmetric fields, such as the Tomimatsu-Sato solutions, are now found 
to be inconsistent with the boundary conditions for colliding plane waves. 

Hassan, Feinstein and Manko (1990) have adapted a generating tech- 
nique of Gutsunaev and Manko (1988) for stationary axisymmetric space- 
times to derive solutions which describe the interaction following the col- 
lision of plane gravitational waves. This technique generates solutions 
in which the approaching waves have variable polarization from a 'seed' 
solution Vo of the equation (10.13) for the colinear case. A solution of the 
Ernst equation (11.17) is then given by 

^_ y^il - a'^ah)t-\- ia{a-\-h)z - {1 + iaa){l - iah) ('12 56) 

(1 — a?ah)t + ia{a + h)z -|- (1 -l- iaa){l — iah) 

where a is a constant and the functions a(t, z) and 6(t, z) are obtained by 
integrating the equations 

{\oga),t = -^[{tz - l)Vo,t + (1 - z^)V,,,] 
t — z 

(log a),, = -^[(1 - t')Vo,t + {tz - l)VoA 
t — z 



(log 6),, = -^[{tz + l)Vo,t + (1 - z^)Vo,. 
t + z 

(log 6),, = ^[(1 - t^)Vo,t + {tz + 1)K,. 
t + z 

The boundary conditions then have to be applied. 



(12.57) 



^ For recent reviews see Kramer (1987), Dodd (1990) and the conference 
proceedings edited by Hoenselaers and Dietz (1984). 



116 Solution- generating techniques 

One particularly important technique that has been developed in 
recent years is the inverse scattering method of Belinskii and Zakharov 
(1978, 1979) which had previously been developed as a soliton technique. 
This method has been applied by Carr and Verdaguer (1983) to plane- 
symmetric space-times, and their results have been interpreted in terms 
of cosmological gravitational waves. The relation between this and other 
techniques has been analysed by Letelier (1987). 

The first solution for colliding plane gravitational waves that was 
obtained using the inverse scattering method was that of Ferrari and 
Ibafiez (19876), who obtained the solution described in Section 10.4 us- 
ing the method in the form developed by Carr and Verdaguer (1983). In 
this paper they took as seed metric the Kasner solution (10.26) which is 
equivalent to the Stoyanov solution (10.29). It may therefore be noted 
that, although the new solution satisfies the boundary conditions for col- 
liding plane waves, the seed solution in this case does not. Non-aligned 
solutions have also been obtained using this method by Ferrari, Ibanez 
and Bruni (19876). 

Another solution-generating technique that has been developed in re- 
cent years is that associated with Backlund transformations. These were 
first applied to general relativity independently by Harrison (1978) and 
Neugebauer (1979). Backlund transformations can also be described as 
double-Harrison transformations, or as quadruple-Neugebaucr transfor- 
mations. The explicit relation between this technique and other methods 
has been discussed by Cosgrove (1980, 1982). 

The effect of a double-Harrison transformation can be deduced by 
any number of competing techniques. Hauser and Ernst (1979, 1980) 
developed a method which involves solving a homogeneous Hilbert prob- 
lem. Cosgrove (1981) has shown how Backlund transformations can be 
conveniently expressed in this formalism. The method has been applied 
to colliding wave problems by Ernst, Garcia-Diaz and Hauser (1988) (see 
also Ernst, 1988), and a class of new colliding wave solutions has been 
obtained. 

It has been observed by Witten (1979) that Einstein's vacuum equa- 
tions for stationary, axisymmetric space-times are equivalent to a form of 
the self-dual Yang-Mills equations. He therefore suggested that it should 
be possible to generate solutions of Einstein's equations by using twistor 
methods to construct self-dual Yang-Mills fields with appropriate sym- 
metries. Such an approach to stationary axisymmetric space-times has 
subsequently been developed by Ward (1983), Woodhouse (1987) and 
also by Woodhouse and Mason (1988), who showed how this approach 
relates to the other solution-generating techniques mentioned above. It 
follows that this approach may also be applied for other types of space- 
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times which contain two commuting KiUing vectors. Woodhouse (1989) 
has appUed these techniques to generate cyUndrical gravitational waves. 
Twistor techniques may also be developed to generate colliding plane wave 
solutions. 

Any of the above techniques may be used to generate new solutions 
of Ernst's equation. Most have been applied mainly to stationary axisym- 
metric space-times. They may, however, also be used to generate solutions 
of Ernst's equation that can be considered as possible to colliding plane 
wave solutions. The difficulty is then to find under what circumstances 
the appropriate boundary conditions are satisfied. Whether or not a new 
solution in region IV satisfies these conditions will generally depend both 
on the particular technique that is employed and also on the initial seed 
solution. 



